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We consider model A dynamics for a quench from the disordered into the ordered phase of 
SU(3) lattice gauge theory and the analogue 3d 3-state Potts model. For the gauge model 
this corresponds to a rapid heating from the confined to the deconfined phase. The expo- 
nential growth factors of low-lying structure function modes are numerically calculated. 
The linear theory of spinodal decomposition is used to determine the critical modes. This 
allows for the Debye screening mass estimation in an effective phenomenological model. 
The quench leads to competing vacuum domains, which make the equilibration of the 
QCD vacuum after the heating non-trivial. The influence of such domains on the gluonic 
energy density is studied. 

Keywords: Deconfining phase transition; lattice gauge theory; spin systems; dynamical 
evolution, spinodal decomposition. 

1. Introduction 

In Ref . 01 121 it is emphasized that in studies of the QCD deconfining transition (or 
cross-over) by means of heavy ion experiments, one ought to be concerned about 
non-equilibrium effects after a rapid heating of the system. This is because the 
heating may be modelled as a rapid quench and the QCD high temperature vacuum 
is characterized by ordered Polyakov loops which are similar to the low temperature 
phase of analogue spin models. The subsequent evolution leads to vacuum domains 
of distinct Z3 triality, and one ought to be concerned about non-equilibrium effects. 

Here we extend the investigation to SU(3) lattice gauge theory. We report pre- 
liminary results about the influence of such domains on the gluonic energy density 
and pressure of pure SU(3) lattice gauge theory. 

The Markov chain Monte Carlo (MC) process provides model A (Glauber) dy- 
namics in the classification of Ref. EH As a time step a sweep of systematic updating 
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with the Cabibbo-Marinari^ heat-bath algorithm and its improvements of Rcf. 5 , 
Elis used (no over-relaxation, to stay in the universality class of Glauber dynamics). 
Although this is certainly not the physical dynamics of QCD, in the present state 
of affairs it appears important to collect qualitative ideas about eventual dynamical 
effects. For this purpose the investigation of any dynamics, which actually allows 
for its study ought to be useful. 

2. Preliminaries 

2.1. Two-Point Correlation Function 

Consider two-point correlation functions defined by 

(«o(0)t4(j))z = TeY. u ^ u ^ + ^i W 

^ a 

l 

where uq(i) is the relevant fluctuation about some average and i denotes spatial 
coordinates. Periodic boundary conditions are used and the subscript L on the 
left-hand side reminds us that the average is taken over the spatial lattice. For 
gauge systems we deal with fluctuations of the Polyakov loop, and for analogue 
spin systems with fluctuations of the magnetization. 

The finite volume continuum limit of Q is achieved by lattice spacing a — > 0, 
N a — * oo with the physical length of the box L = aN a = const. This means that 

(«o(0)«S(j)>i = -3L E a3u °$ u o(^+/) (2) 

a (7 

I 

transforms into 

(u(0)u\R)) L = -L J d\u{r)u\r + R), (3) 
with f—ai, R — aj, u(r) = Uo(i), and so on. 



2.2. Structure Functions 

We define structure function as Fourier transform of the two-point correlation func- 
tion 10): 

F(p) = [ (u(0)u\R)) L d 3 R. (4) 



Occasionally it is also convenient to use structure functions normalized by the vol- 
ume: 

S(p) = -^J (u(0)uHR)) L d 3 R. (5) 
Periodic boundary conditions imply: 

p=—n, (6) 
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where n is an integer vector (0, 0, 0), (0, 0, 1), etc. 
The discretized version of is 



(7) 



Using definition Q and shifting summation on j one arrives (after simple algebra) 
at the expression 



m = h 



Rewriting the product in the exponent as 

2vr - 



api = a — ni — a 

L aN, 



2ir ^ 2?r ^ 
ni = ni. 



and using definition J5J we shape S{p) into the form used in our simulations: 

2 



S(p) 



1 ^— \ I Z7I l 



2m 



N„ 



ni >Uo{i , 



(8) 



(9) 



(10) 



As we let system evolve with time quantity u (i) becomes time-dependent: u (i,t). 
The time t corresponds to the dynamical process, i.e., in our case the Markov chain 
MC time. We consider an ensemble of systems (replica) and dynamical observables 
are calculated as ensemble averages denoted as (...). Then time-dependent structure 
functions averaged over replicas are: 

2 



S(p,t) = (S(p,t)) = 



1 v-^ f 2m _-A 



u (i,t) 



Similarly, 



F{p,t) = L 3 S(p,t). 



(11) 



(12) 



Early time evolution of the structure functions after the quench is governed by the 
linear theory. The linear approximation results in the following equation for the 
structure function: 



dS(p,t) 
dt 



= 2oj(p)S(p,t), 



(13) 



with the solution 

S(p, t) = S(p, t = 0) exp (2w(p)t) , lj{p) > for \p\>p c , (14) 



where p c > is a critical momentum. Originally the linear theory was developed 
for model B t3Hl Details for model A can be found in Ref. |3 
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Fig. 1. The first structure function mode for the 3D 3-state Potts model on Nj* lattices. 



During our simulations the structure functions are averaged over rotationally 
equivalent momenta and the notation S ni is used for the structure function at 
momentum 

2tt 

p= — n where \ft\ = rij . (15) 
L 

The S ni are called structure function modes or structure factors. We recorded 
the modes (including the permutations) n\: (1,0,0), n^'. (1,1,0), n^: (1,1,1), 714: 
(2, 0, 0), n 5 : (2, 1, 0), n 6 : (2, 1, 1), n 7 : (2, 2, 0), n 8 : (2, 2, 1) and (3, 0, 0), n 9 : (3, 1, 0), 
ni : (3, 1, 1), mi: (2, 2, 2), n u : (3, 2, 0), n 13 : (3, 2, 1), n M : (3, 2, 2), n 15 : (3, 3, 0), n 16 : 
(3,3,1), ni7: (3,3,2), ms : (3,3,3). Note that there is an accidental degeneracy in 
length for n§. 



3. Numerical Results 

3.1. Structure Function Modes 

In Figs, m El the time evolution of the first structure function mode S ni after 
a heating quench is depicted (from (3 = 0.2 — > 0.3 in the 3D 3-state Potts model, 
from j3 = 2.0 -> 3.0 in pure SU(2) and from /3 = 5.5 -> 5.92 in pure SU(3) lattice 
gauge theory). Notable in these figures is the strong increase of the maxima S , ™ ax 
with lattice size. In our normalization of S ni non-critical behavior corresponds to 
a fall-off ~ and a second order phase transition to a slower fall-off ~ 

with < x < 3. As the N a — > oo limit is bounded by a constant, our figures show 
that with our lattice sizes the asymptotic behavior may have been reached for the 
3D 3-state Potts model and for SU(2), but not yet for SU(3) lattice gauge theory. 

Our Potts results are from Ref.El Our results for SU(2) are averages over 10,000 
repetitions for the 4 x 16 3 lattice, 4,000 for 4 x 32 3 , 800 for 4 x 48 3 , 340 for 4 x 64 3 
and 106 for the 4 x 80 3 lattice. For SU(3) we rely on 10,000 repetitions of the quench 
for the 4 x 16 3 lattice, 4,000 for 4 x 32 3 and 170 for the 4 x 64 3 lattice. 
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Fig. 2. The first structure function mode for pure SU(2) lattice gauge theory on 4 X iVj lattices. 
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Fig. 3. The first structure function mode for pure SU(3) lattice gauge theory on 4 X lattices. 



The results^of w(fc) fits from a (5 = 0.2 — > 0.3 quench at zero external magnetic 
field in the 3D 3-state Potts model are compiled in Fig. 0] versus n 2 /N 2 . Approx- 
imately, we find straight lines u(k) = ao + ai n 2 with a negative slope a\ and we 
determine the critical momentum k c = 2Trn c /N a as the value where u>(k) changes 
its sign. We find k c w 0.349 from our combined data taken on the 40 3 , 80 3 and 100 3 
lattices^. 

We use the same techniques to determine the critical mode p c in pure SU(3) 
lattice gauge theory. The results are shown in Fig. and indicate ap c = k c = 
2irn c /N a w 0.34 (from the fi gure n 2 /N 2 w 0.003). Relying on a phenomenological 
analysis by Miller and Ogilvie'^"', p c is related by mry — \Z3p c to the Debye screen- 
ing mass at the final temperature Tf after the quench. The relation k c /Tj = N T a k c 
determines Tf and allows us to convert niD to physical units. For our quench we 
have Tf/T c = 1.57 and get 

niD = V3 N T a k c Tf = 3.7 T c . (16) 
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Fig. 4. Determination of k c for the 3D 3-state Potts model on lattices. 







16 ' — i— 
32 

64 > «- 




1 f I 

























0.002 0.004 0.006 0.008 

n 2 /N 2 



Fig. 5. Determination of k c for the pure SU(3) lattice gauge theory on 4 X JVj lattices. 

For pure SU(3) lattice gauge theory T c = 265 (1) MeV holds, assuming a = 420 MeV 
for the string tension, while for QCD the cross-over temperature appears to be 
around T c w 165 MeV, see Ref. 1131 for a recent review. 

In Fig. |31 we observe that not only the heights of the peaks increases with the 
spatial volume, but also the time i max , S™ K = S ni (t max ), which it takes to reach 
them. Whereas S™ x has finally to approach a constant value, i max is expected to 
diverge with lattice size due to the competition of vacuum domain of distinct Z$ 
triality. 

For the Potts models the Fortuin-Kasteleyn (FK) cluster definition can be used 
to exactly remap the phase transition into a percolation model. In case of the 3D 
3-state Potts model the states substitute for the Z3 trialities of SU(3) lattice gauge 
theory. In the cluster language the competition of distinct vacuum domains can be 
made visibleEl. In Fig.[S]we compare the evolution of geometrical and FK clusters for 
a quench of the 3D 3-state Potts model from its disordered into its ordered phase. 
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Fig. 6. Largest geometrical and FK clusters for the 3D 3-state Potts model quenched from j3 = 0.2 
to /3 f = 0.3 on a 40 3 lattice. 



We plot the evolution of the largest clusters for the three Potts magnetizations 
in zero external magnetic field h. While the system grows competing FK clusters 
of each magnetization before one becomes dominant, geometrical clusters do not 
compete. This picture is unfavorable for the use of geometrical clusters of Polyakov 
loops in gauge theories, for which the FK definition does not exist. 

The process of competitions between the largest FK clusters of different magne- 
tization leads for the proper transition (h = 0) to a divergence of the equilibration 
time in the limit of infinite systems, an effect known in condensed matter physics^. 
Potts studies with an external magnetic field show that a major slowing down 
effect survives when h ^ is sufficiently small. As the influence of an external mag- 
netic field on the Potts model is similar to that of quarks on SU(3) gauge theory 
this indicates that the effect may be of relevance for QCD studies of the crossover 
region. 



3.2. Gluonic Energy Density 



Although a satisfactory cluster definition does not exist for gauge theories, the 
underlying mechanism of competing vacuum domains is expected to be similar as 
in the spin models. To study its influence on the gluonic energy e and pressure p 

densities, we calculate these quantities at times t < £ m ax- 

The equilibrium procedure is summarized in Ref.lS HTD fin earlier work EllSI t he 
pressure exhibited a non-physical behavior after the deconfining transition and the 
energy density approached the ideal gas limit too quickly because the anisotropy 
coefficients were calculated perturbatively) . We denote expectation values of space- 
like plaquettes by P a and those involving one time link by P T . The energy density 
and pressure can then be cast into the form 

.,2 



(e+p)/T* = W c N}g- 



2 



r (o) - c T (a)] 



{Pa ~ Pr 



(17) 
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Fig. 7. SU(3) gluonic energy density: (a) with Pq calculated from the time series after the quench 
and (b) using equilibrium values for Pq. 



where Pq is the plaquette expectation value on a symmetric (T — 0) lattice, and 
the anisotropy coefficients ev iT (a) are defined by: 



They are related to the QCD /3-function and can be calculated using Pade fits of 
Ref. As was shown in Ref. ^] this procedure can be successfully used for the 
quench. To normalize to zero temperature, plaquette values from the symmetric 
N T = N a lattice are needed in Eq. (|18|) . As one stays within the confined phase 
on the symmetric lattice its equilibration after the quench is fast. Therefore it is 
enough to use equilibrium values of Pq (at final (if = 5.92) after the quench. This 
is illustrated in Fig. [7| 

In Fig. IHlwe report how the time needed for equilibration increases with lattice 
size. We compare the evolution of the gluonic energy density on lattices 4 x 16 3 , 
4 x 32 3 and 4 x 64 3 . A slowing down of the equilibration is found, which appears 
to be related to the divergence of the time needed to reach the structure factor 
maxima. 

Finally, we compare in Fig. for the N a — 16 lattice the gluonic energy distri- 
bution in equilibrium at (3 — 5.92 with the one obtained after 148 time steps. We 
find a shift towards lower gluonic energies and the width of the probability density 
is slightly broader for the time evolution after the quench than in equilibrium. One 
also has to take into account that the geometry of relativistic heavy ion experiments 
experiments is reasonably approximated by N T /N a — const, N a — > oo, rather than 
by N T — const, N a — > oo. 



and 



(e - 3p)/T 4 = 12N C N$ Ma) - c r (a)] [2P - (P a + P T )\ 



(18) 




(19) 
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Fig. 8. The time evolution of SU(3) gluonic energy density on 4 X N% lattices. 
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Fig. 9. SU(3) gluonic energy density P(e) histograms on a 4 X 16 3 lattice: (a) with competing 
vacuum domains present and (b) after reaching equilibrium. 

4. Summary and Conclusions 

Using the evolution of structure function modes, we identify spinodal decomposition 
as the transition scenario for quench from the disordered into the ordered phase of 
SU(3) lattice gauge theory and the analogue Potts model. We observe an early 
time development of structure factors S n (t) , which is in over-all agreement with the 
exponential growth predicted by the linear theory of spinodal decomposition for 
\p\ < p c . F rom our data the critical mode p c is estimated. Using phenomenological 
arguments Ell p c [ s used to determine the Debye screening mass mo at the final 
temperature T. 

With increasing lattice size N a the time to reach the structure factor maxima 
diverges. Relying on a study of Fortuin-Kasteleyn clusters in Potts models El we 
assume that the reason for the slowing down is due to competing vacuum domains. 
For SU(3) gauge theory these are domains of distinct Z3 trialities. They may be the 
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relevant vacuum configurations after the heating quench in relativistic heavy ion 
collision experiments. We have initiated a study of the gluonic energy and pressure 
densities on such configurations. Our data are consistent with a divergence of the 
equilibration time of the gluonic energy density, similarly to the one observed for 
the structure factor maxima. 

All our results rely on using a dissipative, non-relativistic time evolution, be- 
lieved to be in the Glauber universality class. The hope is that the thus created 
non-equilibrium configurations may exhibit some features, which are in any dy- 
namics typical for the state of the system after the quench. This hope could get 
more credible by studying a Minkowskian time evolution of Polyakov loops and 
finding similar features. Such a study appears to be possible within a relativistic 
Polyakov loop model which was introduced by Pisarski ^1 
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